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Introduction 


Sine ' approximately 1968 there has been an explosive development of 
seconu-order modeling of turbulence (also called invariant modeling or one** 
point closure). This is described in detail in Lumley (1978) » and in sever" 
al other survey articles. Briefly, equations are carried for the various 
second order quantities (the variances and the fluxes), and the unknown 
terms appearing in these equations are modeled as functions of the second 
order quantities, so as to obtain a closed system. 

The approach has been remarkably successful. This is probably because 
more of Che physical mechanisms are carried exactly (unmodeled) in the 
second order equations. Of course, constants are included «diich are adjust" 
ed by calibration of the model against benchmark situations. If the second- 
order models are applied in a situation in which a first-order (mixing 
length) model is adequate, the second-order model usually is not an improve- 
ment. Probably here the additional mechanisms carried exactly are of minor 
importance (or the first-order model would not work), and the crudeness of 
the various closures adopted probably negates the slight additional accura- 
cy. The second-order models will, however, deal with many situations in 
which the first-order models do not work at all (because of the presence of 
the additional mechanisms). 

There are still many ways in which these closures are not adequate, how- 
ever. The constants are by no means universal. In some situations it is 
clear that a physical mechanism has been omitted - for example, two- 
dimensional and axisymmetric jet flows cannot be predicted with the same 


constants (Launder & Mor , 1979). The effect of swirl has not been includ" 
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ed (Lumley, 1980. The closures generally used do not guarantee realLsabil** 
ity: that is, they do not guarantee that quantities that should be non-nega^ 
tive will remain so, or that correlation coefficients will remain less than 
unity in absolute value. 

We feel chat some of these problems can be alleviated by a more formal 
approach to the development of these closures. Although certain of the con- 
cepts of Rational Mechanics have been used from the beginning, they have not 
been applied as extensively or consistently as they could have been. Real- 
izability has not been extensively used to develop closure forms. Many 
closures are essentially ad hoc, and are not tied to a physical model. If 
closures were developed from specific physical models, it should be possible 
to obtain the values of the constants from these models, rather than cali- 
brating Che models against benchmark flows. 

Lumley (1978) attempted Co implement these ideas. In many ways the 
closure presented there is not markedly different from those in common use. 
Some of Che forms were minor improvements; in other cases it was clear what 
was wrung with the forms in common use, but it was not clear how to devise 
an improvement. In one respect, however, the closure presented in Lumley 
(1978) was markedly different: the third moments (the fluxes of the vari- 
ances and fluxes) were obtained from an orderly perturbation procedure about 
a Gaussian equilibrium state. The third moment equations thus obtained were 
minor variations on those suggested on an ad hoc basis by Hanjalic & Launder 
(1972), with one important difference; no addjr - ional adjustable constants 
were introduced. Thus, the transport of fluxes and variances cannot be 


separately adjusted in a calculation. 


Briefly » « turbulence was envisioned whose energy containing scales 
would be Gaussian in the absence of inhomogeneity, gravity, etc. An equation 
was constructed for a function equivalent to the probability density, the 
second moment of vdiich corresponded to the accepted modeled form of the 
Reynolds stress equation. The third moment equations (Stained from this 
(which were thus guaranteed to be consistent) were simplified by the assump- 
tion of weak inhomogeneity. 

The purpose of this paper is to present calculations with this model, 
and interpretations of the results. In the following section we give an 
outline of Che model. The details of the development are given in Lumley (- 
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2. The Mode l 


The dynamic equation for Reynolds stress in an isothermal turbulent flow 
can be written 


*/ * Fi 


( 1 ) 


which states that the substantial change of ’uX“j during convection is 
equal to the production of Reynolds stress from the mean flow plus the 
divergence of the flux of Reynolds Stress and pressure due to turbulent 
fluctuations plus the return toward isotropy of the Reynolds stress tensor, 
olus Che "rapid” chin.je of plus Che (isotropic) dissipation of 

u£uj into heat, respectively. The quantity ^ is the total dissipation. 
At large Reynolds numbers turbulence at the scales at which dissipation 
occurs becomes isotropic, hence, * ^ ^ 

in equation 1. The pressure-strain correlation has been divided into two 
parts, associated respectively with the two Poisson equations 
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2 . d/» . 


( 2 ) 

(3) 


where the first is related to the mean velocity gradient and is linear in 
the fluctuating velocity, while the second is quadratic in the fluctuating 
components. The second is responsible for the return of anisotropic 
turbulence toward isotropy. 

The quantity in the square brackets in equation (1) is a symmetric 
tensor with zero trace, and vanishes if the turbulence is isotropic. Since 
this term acts to interchange energy anong the components vdieo the turbu- 
lence is anisotropic it is natural to express it in terms of the anisotropy 
tensor of the Reynolds stress 
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( 4 ) 


giving 




(5) 


where the inclusion of £ makes the expression dimensionally correct. The 
coefficient can depend on the Reynolds number and the invariants of the 
anisotropy tensor, since must itself be an invariant. Hence. 


C,C3 ,M,Rt) ( 6 ) 

vhece g : ' ^ , M' f ^ 

•nd s . If any one of the velocity components 

vanishes or if Schwarz's inequality between any two components is on the 
verge of being violated, it can be shown that C^-^Z. Hence, to ensure 
realizability, that non-negative quantities, such as the component 
intensities, remain positive and Schwarz's inequality is satisfied we must 
have C ^ > 2. Furthermore it can be shown that 


^ J4T ♦ ^ > 0 (7) 

and is zero only when one eigenvalue of TT^uJ vanishes. Both the 
vanishing of a component and the violation of Scharz's inequality are 
equivalent to the vanishing of an eigenvalue. These considerat ions suggest 
that C} is of the form 

C -. X* < 8 ) 


The function F must reflect the trend that C\*l as Re-^0 and C^-e-Z as 
Resedas indicated by the data of Comte-Bellot S Corrsin (1966). K function 
that satifies these conditions and fits the data of Uberoi (see Lumley & 
Newman 1977) in the raid range of anisotropy is 
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r r <x^/-7.77/?/7/mAA <») 


The condition of realizability is then satisfied with the form of Cj given 
by equations (8) and (9) . 

The rapid term f ^ therefore, 
redistributes energy among the component intensities. Using the Fourier 
transform solution of equation (2), assuming a homogeneous mean field, the 
rapid term becomes 




where 


( 10 ) 


and Sqj is the spectrum of the Reynolds stress. From equation (10) it 
is apparent that the rapid term arises from the interaction of the turbu- 
lence with mean flow velocity gradients. The fourth order tensor must 
sat isfy 



for symmetry, lpiij*0 for incompressibility, 
and when the turbulence is isotropic 






which can be determined directly from equation (11) since the form of the 
spectrum for isotropic turbulence is known. In general this tensor would at 
least depend on the anisotropy tensor and the Reynolds number. We take the 
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tensor to be linear in the anisotropy tisnsor. A form which satisfies all 
the above requirements is 

* ^‘■'% *■ "^v ^ ^ 

There is no information on the variation with Reynolds number and we will 
take C to be a fixed constant. 

The constant C has somewhat different values when evaluated for differ" 
ent flows. Reynolds (1976) found C ■ -0.1 works better for the experiment 
of Tucker & Reynolds (1968) and C * -0.2 works better for the flow of 
Ch«Bpagne« Harris & Corrsin (1978). Reynolds recommends C * -0.15. 
Launder^ Reece & Rodi (1975) use a value of -0.145, which they base on 
homogeneous experiments. Lumley (1978b) suggests C ■ -0.166 to give a 
reasonable agreement with experimental data. In the present calculations we 
use C *■ -0.15. He note that increasing or decreasing C by 10 percent 

produced only slight variations in the results for the wake: if, by changing 
C, better results for the component energies were obtained, then poorer 
results in comparison with experimental data were obtained for the shear 
stress . 

By formulating a dynamic equation for the cumulant C^jj^ ^corre- 
sponding to the triple velocity correlat ion^ and performing an order->of-mag- 

nitude analysis for the case of weak inhomogeneity (homogeneous turbulence 
is observed to be approximately Gaussian in the energy containing eddies and 
departure from Gaussian behavior is associated with inhomogeneity, that is, 
non-zero values of U£UjU^ are fluxes and are necessarily non-Gaus- 

sian) it can be shown (see Lumley 1978a) 


s - 

^ j 'll. 






j I 





^ 

f to 



( 14 ) 
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where 






The quantity C^ is the same coefficient given by equation (8) for the 
return to isotropy. 

For the pressure transport, (puj,£ ^ pu£,j)/^, we begin from 

equation (3) for the second (return to isotropy) part of the pressure for a 
hmogeneous flow and obtain by use of Fourier transfoms 




(15) 


where S^jj^ is Che spectrum of u^ujuj^. 


Defining 




(16) 


we have *)^»®try, • 0 from 

incompressibility, and “ *V*^>*V*’ most general Linear form for 

contains five coefficients; however, applying the various conditions, 
we can determine them all: 




then 




(18) 


where qZu^ is given by equation (14). If UfeO^ the expression vanish- 
es, and realizability is satisfied. 

To complete the set of equations Che dissipation must be calculated. 
For this we use a convection-transport equation with a production 


/destrucciun term 


where a dimcneionlese invariant function. It it reatonable to ttsuM 
that depends on the Reynolds stress, the mean velocity gradient, the dis~ 
sipation and the viscosity. Hence, 

This nnist now be a function of the invariants that can be constructed from 
these quantities. We suke the stipulation that the mean velocity gradient 
does not appear without the anisotropy, because we do not expect to change 
the level of the dissipation by a change in the mean gradient if the turbu* 
lence is isotropic. There are a large niimher of invariants that can oe 
formed between the mean velocity gradient and the anisotrppy and we only 
give a few*. 


k _ 




( 21 ) 


Assuming the velocity gradient to be small we can expand this function in a 
power series 


^ /f *# ^y/ f f ••• 


( 22 ) 


where the coefficients are functions of the invariants of the anisotropy 
tensor and Reynolds numbers. 

The coefficients in the above expression must be determined from experi*' 
mental data. Luraley and Newman (1977) found that * 14/5 for two 
limiting cates with no production from mean velocity gradients: the final 
period of decay for small Reynolds numbers and for one-dimensional 
turbulence. With these considerations the following formula approximates 
the data of ComteBellot & Currsin (1966) 
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t I*l/f t- 4Hft[ -i.ti (23) 

For the production term we uee Che value ■ 2..0 ae given by Reynolda at 
was determined from the nearly homogeneous shear i!Lows of Tucker & Reynolds 
(1968) and Chaispagne, Harris & Corrsin (1966). Assuming the anisotropy to 
be not too large we ignore they^-term in equation (22). 

To determine the transport flux €4f^ for dissipation we note Chat the 


eddy viscosity not vary greatly across sxich of Che width of free 

shear flows. Also ^^)/T must remain finite as both ^ and q2 vanish ( 
at the edge of a wake or jet.) Presuming that /^^const ant we have 


r ^ f 9 % 




(2A) 


Drawing an analogy from this relation we can write for the transport 




(25) 


where we have included pressure transport in the total transport. This 
relation stipulates that as the energy transport vanishes so must Che dissi~ 
pat ion transport. Using equation (18) for the pressure transport we c^tain 


Us 


nr- 

ing equation (14) for but replacing 


(26) 


we uhtain for the dissipation transport flux 


( 27 ) 
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Although obtained on aowcuhat tenuoua grounda this raault givaa about the 
right iMgnitude for the diaaipation tranaport flux and it haa the advantage 
that no neu conatanta are introduced. 

2. Calculation o f ^ he Wake 

The iaothemal far'*wake downatreae of an <d>ject In an initially unifona 
flou providea a relatively aiaple yet comprehenaive teating ground for a 
turbulence model. The aiaplif ication ariaea froai the fact that aome dia~ 
tance dounatre«a of the object the velocity defect ia email and the convec' 
tive term can linearised. However, all poaaible ingredienta for an iaother* 
mal turbulent flow are preaent; convection, mean velocity gradienta, aniaot* 
ropy, production, and redistribution due to rapid terms and tranaport. 

In terms of the dimenaioneaa velocity defect U ’“^iV^^nd dimen* 
sionl*»^3 citordinaCi** x y ■ X2/I where d ia the mumentue thickness 
of the wake, the component equations for the two*diaMnaional wake are given 
in Appendix 1 and those for the axiaymmetric wake in Appendix 2. The only 
mean flow component of significance is U and except for the eonvfM:tion terms 
all gradients in the flow direction can be neglected. 

Since the equations are parabolic in the x-direction and boundary valu^ 
in the y-direction we can solve them by marching in the x-direction from an 
initial station, solving a boundary value problem in the y**direction at each 
x-location. A simple implicit scheme was formulated, using a simple back- 
ward difference in the x-direction with central differences in the y-direc- 
tion. At each x*position the cquati<ma 'or the teynolds stresses are solved 
sequentially. The diffusion term in the variable in question, along with 
part of the return-to-isotropy a^ convection terns, are treated implicitly. 
The other diffusion terms, production, rapid terms, etc. are treated 
explicitly and are evaluated at the upstream station. Finite differencing 
the implicit terms gives a sec of algebraic equations whose coefficient 
matrix is tridiagonal and which is readily solved by an elimination 
procedure. After the Reynolds stresses are computed the dissipation 
equation is solved in a similar manner. Then the mean velocity defect U is 
determined from the momentum equation. 

The boundary conditions at the centerline are syasaetry conditions on all 
the variables except for uv whose value is zero there. At the outer edge 
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all variables approach zero. However, since some terms in the equations 

contain or ? in the denominator, it is necessary to add a small nuodier, 

taken as a small fraction of the center values, to q^ and ^ in order to 
avoid dividing by zero in the calculations. The outer edge was taken at 
yj^ » 4it where X i* the position where 0/U^“ exp (”.5). For the first sta- 
tion 50 intervals were taken across half the wake. As the wake grows more 

intervals are added corresponding to ® 4^. until there are 100, at tdiich 
point y is doubled and every other point is retained so that there are again 
50 intervals. 

Since we are interested in obtaining the far field self-preserving wake, 
the initial conditions specified at x ®* 1, were chosen soaiewhat arbitrarily. 
The similarity solution with constant eddy viscosity given in Tennekes & 
Luffiley ^1972) was used for the velocity profile, ” exp (-0.5 t 

where s i.M'H and ( ' 1?7 for a circular cylin- 
der and \Ja s 0.*?^ and 1, « for a sphere. The key 

input parameter is the value of the momentum integral which is J Udy " 1 

.oo . 

for two dimensions and J Urdr =* 1/8 for the axisymmetric wake. The 
values of these integrals were continuously monitored in the calculations to 
ensure their constancy. The initial shear stress profile uv was also deter- 
mined from the eddy viscosity solution. The turbulent velocity is given by 

3 -7-^ uv and the dissipation is roughly the production for y>Jt, 
Tennekes & Lumley (1972). For 0 < y^jP, q^. and € are taken equal to the 
values at y X • “nie component energies were taken to be one third of the 
total turbulent energy. 

For the two-dimensional case the results are presented for x " 1000, 
after 500 steps, where the profiles are self preserving when scaled by 
and X • Scaled results at x 2000 were less than a few percent different 
than those at x " 1000. Changes in the initial conditions produced little 
difference in the self-preserving behavior of the solutions. Results for 
the axisymmetric case are presented for x ** 500, after 500 steps. A numer- 
ical instability occurred around x ** 700 in attempts to compute to larger 
distances. Even so, the profiles were nearly self-preserving at x *» 500. 
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Results and Discussion 


The calculated predictions for the mean velocities and Reynolds stresses 
are compared with experimental results of Chevray & Kovasznay 1969 (for the 
wake behind a flat plate) and Townsend 1956 (for the wake behind a circular 
cylinder) in Figures 1-5. Overall the agreement between predictions and the 

—V J 

flat-plate experimental results are quite good. The data for and 

for the cylinder are quite different from the predictions. We will pre- 
sent below our -explanation for the discrepancy between the cylinder data and 
the calculated results. We include the cylinder both to contrast the behav- 
ior of wakes behind different objects and to indicate the trends of certain 
turbulent quantities that have only been measured for a cylinder. It is 
seen that the calculated mean velocity profiles are somewhat flatter than 
those measured for the flat plate wake, and the calculated shear profile is 
slightly lower than the measured. For self preservation these quantities 
are directly related by 




(29) 


Probably there is a little too much diffusion in the dynamic equation for 
the Reynolds stress, thus giving too small a peak in the shear profile and a 
flatter velocity profile. We note that Townsend's shear data does not obey 
the self preserving relation (29). 


Figure 6 shows the calculated transport flux profiles 
and as compared to the measurements of Townsend. 


Although, as we 


shall argue below, the predictions should not compare well with the cylinder 
wake, the trends compare reasonably and the magnitudes are about right. 
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T^ble I gives the overall characteristics of the plane wake. It is seen 
that there is good agreement between the predicted wake growth and center- 
line velocity decay and the flat plate wake data. We suggest that this is 
so because the flat plate wake is created rather smoothly from the trailing 
edge boundary layers on the top and bottom of the thin plate, without vortex 
shedding. These conditions seem to us to be a better match to the theory, 
which does not specifically include transport by large eddy structures. 

The cylinder, on the other hand, has alternate vortex shedding. As the 
flow ages these vortices Lose their identity to a certain extent, but a dis- 
tinct large-scale structure remains. It is generally accepted that far 
downstream, where the flow is in dynamic equilibrium, the wake becomes inde- 
pendent of everything except the initial momentum thickness. However, this 
universal equilibrium wake has not been observed experimentally, perhaps 
because measurements have not been carried out to large enough distances. 
What have been observed experimentally (Bevilaqua & Lykoudis 1978) are 
apparently self-preserving wakes, which obey the similarity laws within 
experimental error over the range of observation, but which are different 
for different objects, e.g. , cylinder versus flat plate. This apparently 
self-preserving mode persists with seemingly little change for considerable 
distances, as is evidenced by Townsend's cylinder measurements which exhibit 
the same behavior to x/d > 1000. It seems possible that the effects of the 
near-field large scale structure (characterized by meandering of the entire 
wake, with strong intermit tency, and the occasional appearance of non- 
turbulent flow on the centerline) persist to great distances, the wake 
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continuing to evolve, but so slowly as to be undetectable c»ver the range of 
observation. It is also possible that the flat plate wake will slowly 
evolve its own large scale structure far downstream. This could explain the 
differences between the plate and cylinder wakes: the plate having little or 
no near^field large-scale structure quickly attains an apparent equilibrium 
similarity mode (from which it may be slowly evolving toward a 
state) whereas the cylinder quickly attains a different apparently self- 
preserving mode, %diich may also be slowly evolving toward a universal state; 
the two will not reach their common universal state until far beyond the ex- 
perimental range. 

Lc was stated in the description of the numerical procedures chat Che 
initial conditions on centerline velocity and characteristic wake dimension 
were chosen consistent with the measurements for the cylinder given in Table 
1. Yet shortly into the the calculations, the wake characteristics tended 
toward Che asymp^tic calculated values which are close Co those given for 
the flat plate wake. In order that Che calculated constants given in Table 
1 be determined independent of Xq we used and 0“^^/ 

(k » which corresponds to the relations Xs and 

,\*^espectively. These constants closely assumed their asympototic 
values around x s 100. Furthermore, changes in initial conditions, energy, 
dissipation levels, etc., (except for Che initial momentum thickness) had 
little effect on Che ultimate self-preserving state, affecting only the rate 
at which it was attained. Despite such changes, the results tended to seek 
a unique self-preserving state, which we may speculate is the universal sim- 
ilarity equilibrium state without large scale structure. 

Calculations carried out with the model of Launder, Reece, & Rodi 4975) 
also produced good results for the self-preserving profiles when compared to 
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the experimental results of Chevray & Kovasznay. In their model » 
Cp ■ 0.055 for the diffusion coct.?icient in the dynamic equations for the 
Reynolds stresses, C “ -0.0145, ■ 0.075 for the diffusion function in the 
dissipation equation, ^0 » 3.8 and 1 ■ 2.88. The closure formulations 
are those given by those authors. Ue found good results for the centerline 
velocity variation and the growth rate of the characteristic lateral dimen- 
sion. However, in Hanjalic & Launder (l980^ it is stated that there is no 
explanation why the theory predicts a growth rate vihich is 35Z less than 
measured. These authors have- contrasted their predictions with measurements 
for the circular cylinder. 

Rodi ^1975^ organized data from ten wake studies: six circular cylinders, 
one normal plate, two airfoils and one aligned plate. Although there are 
some differences^ most of these studies show consistent results for the 
growth constant. Rodi points out that the Gievray & Kovasznay flat plate 
experiment was not fully-developed in the distance in \<hich they took meas- 
urements. It is true that the scaled turbulent intensity and shear were 
still changing, although the mean velocity appears to have reached self- 
preservation (the mean velocity field usually reaches self-preservation well 
before the mean turbulence field does) . However, we feel that the wake of 
Chevray & Kosvasznay is closer to some sort of (at least temporary) equilib- 
rium chan it is to the cylinder wake, to which it does not appear to be 
tending. 

The form of the equations for the axisymmetric wake are much the same as 
those for the plane wake except that Che diffusion-transport terms are much 
more complicated. These terms have been worked out and are presented in 
Appendix II. Results for the round wake^ calculated by procedures explained 
for the plane wake, are presented in Figures (7-11). Figure 7 shows excel- 
lent agreement for Che mean velocity between Che calculated results and 
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Che meaauremeaCs of Bevilaqua & Lykoudis (l978), for Che wake downsCreani of a 
porous disk 2.54 cm in diameter (normal Co Che stream) and of Qtevray (1969) 
for Che wake downstream of a 6 x 1 prolate spheroid with a length Reynolds 
number of 2.75 x 10^ (so chat the boundary layer on Che surface of Che 
spheroid was turbulent). Both of these wakes were relatively smooth in Che 
near field with no evidence of vortex shedding or large scale eddy structure 
from flow-visualizacion experiments. Hence, we hypothesize these wakes 
would quickly attain an apparent equilibrium similarity. The agreement 
between the measured and calculated normal Reynolds stress is also fair. 
However, comparisons for and for Chevray*s spheroid 
shown in Figures 9-11 are not so good. The reason for this is probably chat 
the measurements were taken at x/d ■ 18, a position not far enough down- 
stream fur the turbulence quantities to become self-preserving. 

Table 2 gives the overall characteristics for the various wakes. Ic is 
seen that the calculated wake is close to the porous disk in centerline 
velocity decay and growth rate where iMasurements were carried out to a dis- 
tance of x/d “ 100, The calculated growth and centerline velocity decay is 
also not too different from those for the spheroid wake. Note that the 
behavior of the wake behind a sphere is quite different from the others: the 
vortex shedding, meandering of the wake and large scale structure appears to 
leave a significant influence on the downstream wake behavior. The differ- 
ence in wakes and the possible existence of different states of self preser- 
vation is the subject of the paper by Bevilaqua & Lykoudis. We only argue 
here that for proper comparisons between experiment and this theory (without 
large scale structure) the experiments '■usi also lack large-scale struc- 
ture. 

As a closing point we must mentiC'i ihe discrepancy that has existed over 
the years in predicting round and plane jets with turbulence models. A 
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theory with a set of constants that worked well for the plane case did a 
poor job of predicting the round jet. We did not find a large discrepancy 
for the wake when the calculations were referred to the experiments lacking 
in large-scale structure. Possibly our more general formulation of the 
return-to-lsotropy coefficient, the transport terms, the pressure transport, 
etc., has lead to a more nearly universal theory capable of handling plane 
and three-dimensional flows in the absence of large-scale structure. 

For the treatment of flows containing large-s.cale structures (presumably 
the majority of flows important in technology and nature) it will probably 
be necessary to introduce the structure explicitly by some form of stability 
analysis. One of us (JLL) has recently made suggestions along these lines 

h 

(Lumley, 1981^). 
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Appendix 1 

The equations for calculating the plane wake are presented here. 
They are written in terms of the mean flow velocity defect II ■ (Uios *" 
Ui)/Ujo»; the non-dimensional coordinates x ■ xj/;^ and y ■ yi/i^f where xj 
is in the mean flow direction and yj is measured across the shear layer; 
and the non-dimensional turbulent velocity components u ■ u]^/U|os» w “ 
U2/U190, and w - 03/11^00. The usual boundary layer assumptions and wake 
s imp I if ic at ions are employed. 


Momentum equation: 


, ' 3-177 

Reynolds stress equations: 
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Appti^ix II 


The equations in non-dlnens tonal forts for the axiaysMMtric wake are 
presented in teras of the mean velocity defect (defined as for the plane 
wake) and the velocity components u ■ v ■ U2/tI|oe ^ 

for the axial (x - xifp), radial (r - X2/tf*) «nd angular directions 
respectively. 

Momentum equation: 


15 " t * 1 — //-Ml/*) 
rir< 


Reynolds stress equations: 


If - ‘ "-'I'M/"- ‘ 
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Dissipation: 




— -r«. 


f 



t 



, ^ mV* ‘b'O' 

V', ^ ^ jr 


Chevray & 
Kovasznay 
(flat plate) 


(i/ui)* 

((lU-U(x,0)) 

/uj- 

l(x-xo)Ml^/2 

Rf 

0.277 

1.44 

1 

10.4 

0.216 

2.06 

19.1 


Calculated 


0.201 


2.10 


20.9 










ROUND WAKE CHARACTERISTICS 



Chevray 

(spheroid) 


Bevilaqua & 
Lykoudis 
(Sphere) 


Bevilaqua & 
Lykoudis 
(Porous Disk) 



urn- 

l(x-xo)Mll'3 

[(lU-U(x,0)) 

/lUl* 

l(x-xo)/«2/3 

RX 

0.272 

2.09 

23.05 

0.452 

0.575 

3.82 

0.246 

2.31 

28.17 

0.220 

2.48 

33.8 













Figure Captions 


Figure 

Figure 

Figure 

Figure 

Figure 

Figure 

Figure 

Figure 

Figure 

Figure 

Figure 


1. Notmelized mean velocity, plane wake. 

2. Normalized streamwise energy, plane wake. 

3. Normalized energy out of plane of wake. 

4. Normalized cross'-stream energy in plane of wake. 

5. Normalized shear stress, plane wake. 

6. Normalized triple correlations, plane wake. 

7. Normalized mean velocity, round wake. 

8. Normalized axial energy, round wake. 

9. Normalized radial energy, round wake. 

10 Normalized azimuthal energy, round wake. 

0 

11. Normalized shear stress, round wake. 
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